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Abstract 

We study the effects of dark matter annihilation during and after the cosmic 
recombination epoch on the cosmic microwave background anisotropy, taking into 
account the detailed energy deposition of the annihilation products. It is found 
that a fairly stringent constraint on the annihilation cross section is imposed for 
TeV scale dark matter masses from WMAP observations. 



1 Introduction 



Although the existence of dark matter in the Universe is promising, its observational 
evidences rely only on its gravitational properties. Despite tremendous efforts, particle 
physics natures of the dark matter remains unrevealed. 

Many possible ways to detect the dark matter were proposed so far pQ. Direct detec- 
tion experiments search for the dark matter scattering signatures off nucleons. Indirect 
detection experiments search for characteristic cosmic-ray signals originated from dark 
matter annihilation inside and/or outside the Galaxy. Recently the PAMELA satellite 
reported an excess of cosmic-ray positron flux around the energy 10-100 GeV [2], and it 
can be interpreted as a signal of dark matter annihilation. Results of the Fermi satellite [3] 
on the cosmic-ray electron plus positron flux seems to be consistent with the dark matter 
interpretation of the PAMELA excess. If this scenario is correct, however, a fairly large 
annihilation cross section, typically (av) ~ 10 _24 -10 -23 cm 3 s _1 is required, which is two- 
or three-orders of magnitude larger than the canonical value ((crv) ~ 3 x 10~ 26 cm 3 s _1 ) 
reproducing the correct dark matter abundance under the standard thermal freeze out 
scenario. It was pointed out that the existence of dark matter with such a large annihi- 
lation cross section is constrained from many observations : gamma-rays from Galactic 
center [I] or from extragalactic background [5j [6], neutrinos [7], anti-protons [5], syn- 
chrotron radiation, and even from big-bang nucleosynthesis (BBN) [9]. More or less, all 
of these constraints so far seem to be comparable to one another. Some models can be 
killed by some observations, but there is a possible way out. For example, the Galactic 
gamma-ray constraint is often very stringent but sensitive to the density profile of the 
dark matter halo. It can be evaded by assuming the flat cored profile. 

In this paper we propose a more rigorous way to constrain dark matter annihilation 
cross section. Although dark matter annihilates with each other throughout whole history 
of the Universe, let us now focus on its effects around the recombination epoch. Dark 
matter annihilation injects high energy particles into thermal plasma consisting of pho- 
tons, electrons, protons, neutral hydrogens and heliums. High energy particles interact 
with background plasma and some fraction of their energy is used for the ionization of 
neutral hydrogens. Therefore, dark matter annihilation affects the recombination history 
of the hydrogen resulting in higher residual ionization fraction (x e ) than the standard 
value. Even a slight shift of x e can cause a drastic effect on the cosmic microwave back- 
ground (CMB) anisotropy, and hence the annihilation rate is constrained from the precise 
measurements of the CMB anisotropy. 

Similar analyses were carried out in several papers [TUJ [11], [12], [13] , especially in con- 
nection with recent anomalous cosmic-ray positron/electron fluxes [H] |T3], [121 [13 HE]- 
Many of those analyses are based on the study in Ref. [TU] where effects of high-energy 
electron/photon energy injection due to late decaying particles on CMB were investigated 
using the extrapolation of the result of Ref. [20] for the energy fraction going into the 
ionization, heating and excitation of the intergalactic medium, which was valid only for 
low energy electron/photon injection. Ref. [2TJ derived a formalism to extend this to 
high-energy electron/photon energy injection, but there the Hubble expansion effect was 



1 



neglected. Therefore, we further extend the formalism of Ref. [2T] to correctly handle the 
Hubble expansion, and apply it to the study of the CMB constraint on the dark matter 
annihilation. (See Ref. [18j for a different approach to compute the details of the energy 
deposition.) 

This paper is organized as follows. In Sec. [2] we generalize the formulation of Ref. [21] 
to include the effect of Hubble expansion for calculating the energy deposition from high 
energy photons and electrons. In Sec. [3] we apply the results of Sec. H]to derive the effects 
of dark matter annihilation on the CMB anisotropy, and constrain the annihilation cross 
section. In Sec. H] we will mention the effect of dark matter clustering. Sec. [5] is devoted 
to the conclusions. 

2 Energy deposition from dark matter annihilation : 
computational methods 

In this section we describe our method to calculate which fraction of the injected electron 
and photon energy goes into the ionization, heating and excitation of the intergalactic 
medium. (The case of another injected particles, such as a muon, will be discussed in 
the next section.) Our method is closely based on Ref. [21] where the effect of Hubble 
expansion was neglected. We extend the formalism to include the Hubble expansion. 

Let Ei, E%, . . . , En be a discrete set of energies of particles (Ei < Ej for i < j). 
Suppose that an electron in the energy bin % is injected at the redshift z. This electron 
interacts with background particles and redistribute its energy toward lower bins in the 
form of both electrons and photons. Then similar energy degradation processes occur for 
each low energy electron/photon, and again similar processes are repeated until z = 0. 
During this cascade processes, some of the initial injected energy is used for the heating, 
excitation and ionization of the intergalactic medium. In order to deal with this entire 
process, let us also introduce a discrete set of the redshift, Z\, z 2 , ■ ■ ■ , z m (zi < Zj for i < j). 

We define xti^h z m, z n ), Xhi^h z m, z n ) and xli^i, z m , z n ) as fractions of the initial en- 
ergy Ei at z m for an electron (a = e) or photon (a = 7) which goes into ionization, heating 
and excitation at the redshift bin of z n , respectively. We also introduce Xz(Ei, z m, z n) as 
a fraction of the initial energy lost by the Hubble expansion at z n . The remaining en- 
ergy, which is not yet used for the above processes is denoted by Xr{^ii z mi z n)- These 
Xa(Ei,z m ,z n ) (a = i,h,e,z,r) are the quantities which we want to know. In order to 
evaluate them, the following probabilities are necessary. P a (Ei, Ej, z m ) denotes the prob- 
ability for an electron (a = e) or photon (a = 7) with energy Ei is transferred to another 
energy bin Ej by scatterings at z m . The background electron/photons are scattered up 
through the same process, a + 6(BG) ({a, b} = {e, 7}) and the possibility for this is 
represented as P^ G (Ei,Ej,z m ). The possibility for an electron/photon in energy bin % 
redshifts away is denoted by Pjj(Ei,z m ,z m -i). There are two relations which can be 
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easily understood, 

£ P a (E l , Ej, z m ) + P H (Ei, Ej, z m ) = 1, (1) 
P a (Ei, Ej, z rn ) = ^2 Pbg(Ei, Ei ~ Ej, z m ). (2) 

6=e,7 

Let us first consider the case z m = z n . In this case Xa{Ei, z m , z m ) are given by 

Xi,e,h(Ei, z m, z rn) = ^ ] \_P {Ei, Ej , Z m )Xi,e,h(Ej , Z m , Z m ) 



j<i 

b b E 

+Pbg(Ei, Ej, z m )Xi, e ,h(Ej, z m , z m )] -^f , 

Xz(Ei, Z m , Zm) ^ ^ [P {Ei, Ej, Zm)Xz{Ej , Z m , Z m ) 
j<i 

E ■ A E' 

J rPBG\Ei, Ej, Z m )Xz{Ej, Z m , Z m )\ — + P^Ei, Z m , Z m -l) ~ , 

&i &i 

where A Z E { = E^l - (1 + z m -i)/(l + z m )}, and 

Xr{Ei, Z m , z rn) ^ ^ [P (E{, Ej, Z m )Xr{Ej , Z rn , Z rn ) 
j<i 

^PBciEi, Ej, Zm)Xr(Ej, Z m , Z m )] — f + Pjj(Ei, Z m , Z m ~\) — • 



(3) 



(4) 



(5) 



The energy conservation is checked by J2 a Xa{Ei, z m , z m ) = 1. When z n < z m , Xa{Ei, z m , z n ) 
are evaluated by 

Xa(Ei, Z m , Z n ) = ^ ^ [P {Ei, Ej, Zm)Xa(Ej, Z m , Z n ) 
j<i 

+P£ b G (Ei, Ej, z m )x b a (E J , z m , z n )] || (6) 

Ei — A z Ei 



+ Pff(Ei, z m , z m -i)Xa (Ei — A z Ei, z m -i, z n ) 



Ei 



The energy conservation is also satisfied (%2 a Xa(Ei, z m , z n ) = Xr(Ei, z m , z n+ i)). There- 
fore, once we know all Xa(Ej, z\, Zk) with Ej < Ei and zi < z m , we can compute 
Xa(Ei, z m , z n ) using the above equations. 

The remaining task is to determine the transition probabilities : P a (Ei, Ej, z m ), 
Pff(Ei, z m , z m -i) and P£ b G (Ei, Ej, z m ). They are related to the collision frequency of ele- 
mentary processes. We take into account all the relevant interactions. For an electron, 
the excitation and ionization of and collision with the hydrogen atom, Coulomb collision 
with the background electrons, and inverse Compton scattering with CMB photons are 
considered. For a photon, photoionization of the hydrogen atom, Compton scattering 
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with background electrons, pair cration in nuclei, photon-photon scattering and double 
photo pair creation with CMB photons are considered. For simplicity, we neglect 4 He and 
assume that all baryons are H. 

Since the energy loss processes contain both "continuous loss" where the energy loss 
per collision is smaller than the energy bin AEi, and "discrete loss" where the energy 
loss is significant, we need to distinguish them for a numerical purpose. The collision 
frequency for a discrete loss can be given by 



dE, 



31 



(7) 



where n t is the number density of the target particle, v p is the particle velocity, and 
da^iEj/dEj is the differential cross section. The collision frequency for a continuous 
loss is expressed as 



v a c (E h EU> 



AEi 



-dE 
dt 



Similarly, we define the "collision" frequency for the Hubble expansion as 



1 



A z Ei 



-dE 
~dF 



H 



H{z m )Ej 
A-,Ei 



(8) 



(9) 



where H(z) is the Hubble parameter at the redshift z. Using these quantities, transitions 
probabilities are expressed as 

z/ c a (^,E i _ 1 ) + ^(E i ,E i _ 1 ) 



P a {Ei, Ei_ ly z m ) 
P a {E h E h z m ) 

PlI\Eii z mi z m~l) 

where v% = + z/J 7 . We also obtain 
Pbg(Ei, Ei - Ej, z. 



v t{Ei, Ei-i) 




E k ) 


+ v a H {Ei) 










vZ{E h Ei_ x ) 


+ J2 k ^(Ei, 


E k ) 


+ " a H (Ei) 




" a H (Ei) 






K(Ei, Ei-x 


) + E k K(E t 


E ki 


+ v a H (Ei) 



v- c {E h Si_0 + E k "d(Ei, E k ) + iS H (Ei 



(10) 



The elementary processes contributing to P^q are as follows : collisional ionization of 
the hydrogen atom for (a = e,b = e), inverse Compton scattering for (a = e,b = 7), 
photoionization, Compton scattering, photon-matter pair creation, double photon pair 
creation for (a = 7,6 = e), and photon-photon scattering for (a = 7, b = 7). All the 
relevant cross sections and energy loss rates are summarized in Ref. [21] . and we do not 
repeat here. 

The above methods allow us to compute Xa(Ei, z m , z n ). In Fig. [1] the differential 
ionization, heating and excitation fractions dxfl e (E , z' , z) / dz are shown by red-dotted, 
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Figure 1: The differential ionization, heating and excitation fractions dxfl e {E , z' , z) / dz 
are shown by red-dotted, black-solid and dot-dashed-blue lines, respectively. Thick (thin) 
lines correspond to z' = 1000 (100), and we have taken E = 1 TeV. 



black-solid and dot-dashed-blue lines, respectively. Thick (thin) line corresponds to the 
injection redshift z' = 1000 (100), and we have taken E = 1 TeV. Here no reionization from 
astrophysical sources is assumed for simplicity!]] Now we are in a position to accurately 
calculate the effects of dark matter annihilation on the CMB anisotropy, as we will see in 
the next section. 



3 Effects of dark matter annihilation on the CMB 
anisotropy 

Now let us apply the results of previous section to the calculation of CMB anisotropy. 
In order to take into account the effects of extra ionization and heating source from dark 
matter annihilation, we should add the following terms for the evolution of the ionization 
fraction (x e ) and the gas temperature (T&), 



dz 



J DM 



dz' 



rr 



[(z')(av) F m x d X T\E,z',z) 



H(z')(l 



z' 



n H (z') E Ry 



dz 



(11) 



1 In this calculation, the background ionization history is fixed and the extra contribution to the 
ionization is regarded as a small perturbation. We have checked that this approximation does not alter 
the following result at all, since otherwise the resultant cosmology significantly deviates from the standard 
one and already excluded. 
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Figure 2: History of the ionization fraction as a function of redshift. The black solid 
line corresponds to the standard recombination history, without dark matter annihilation 
effects. Also shown are the cases of dark matter with annihilation cross section (crv) e + e - = 
10 -24 and 5 x 10 -24 cm 3 s _1 with the dark matter mass rn x =l TeV. 



where E Ry = 13.6 eV is the Rydberg energy, m x and n x are the mass and number density 
of the dark matter particle, n# is the number density of the hydrogen atom. 



H(z')(l + z') 3n H (z' 



dTb f dz' 

~dz 

U/6 J DM 

Here we have defined 

dxZ\E,z',z) _ 

dz / m v dE dz 



2n\(z'){av) F d X [ F) \E , z' , z) 
■m v — 



d,z 



(12) 



m x 



dN^dxtl(E,z',z) dN^d X ^(E,z' 



(7) W V W/ 



+ 



dE 



dz 



(13) 



where dNp'^ /dE denotes the spectrum of the electron and photon produced per dark 
matter annihilation into the mode F, and (<jv)f denotes the annihilation cross section 
into that mode. We have included these terms in the RECFAST code (25], which is 
implemented in the CAMB code [23J for calculating the CMB anisotropy. Here and 
hereafter, we fix the cosmological parameters to the WMAP five year best fit values |24j . 
The reionization optical depth is also fixed to be the best fit value and need not be 
reevaluated when the dark matter annihilation effect is included, since it depends only 
on the reionization history at low-redshift. It is noted that the energy integral in ffTB"]) for 
given final states F can be performed before solving the evolution equation once we have 



6 



tables of dxfl/ dz and dx\h/dz, as long as the produced particle F decays so quickly as not 
to lose energy by the interaction with surrounding medium, which is a valid assumption 
for any unstable standard model particle. 

Fig. [2] shows the evolution of the ionized fraction with and without the effects of 
dark matter annihilation. The black solid line corresponds to the standard recombination 
history, without dark matter annihilation effects. Also shown are the cases of dark matter 
with annihilation cross section (av) e + e - = 10 -24 and 5 x 10 -24 cm 3 s _1 with m x =l TeV. It 
is clearly seen that energy injection from dark matter annihilation works as an extra source 
of ionization. Then it is not hard to imagine that this modified recombination history 
affects the spectrum of CMB anisotropy, since large x e implies large optical depth. 

Fig. [3] show the multipole coefficient of the CMB anisotropy (Ce) of TT, TE and EE 
modes, respectively. The cases of {av) e + e - = 10 -24 and 5 x 10~ 24 cm 3 s _1 with m x =l TeV 
are also shown. It is seen that as the dark matter annihilation cross section is increased, 
the spectrum is more suppressed and it will deviates from the observed data. 

In Refs. [TBI [TT] , the optical depth was calculated for comparison with the CMB ob- 
servation. In this paper, we have performed a x 2 analysis in order to compare the results 
quantitatively with the real data, using the WMAP likelihood code [25] . In Fig. HI result- 
ing 99% C.L. constraints on the dark matter annihilation cross section are shown. The 
solid (dashed) line corresponds to the case where dark mater annihilates into n + fi~ (e + e~). 
It is seen that CMB measurements gives stringent upper bound on the annihilation rate. 
Note that, however, in this analysis we have fixed cosmological parameters as WMAP 
five year best fit values. It may be expected that there are some degeneracies between 
dark matter annihilation cross section and other cosmological parameters, as shown in 
Ref . [TUl [T5] . Thus the bound will be somewhat loosen when the full Monte Carlo analysis 
is performed. We will give more detailed analysis on this subject elsewhere. 

4 Effects of dark matter clustering 

So far, we have assumed that the dark matter distributes homogeneously in the whole 
Universe. Of course, this is not true especially for the low redshift Universe, since the 
structure formation of the dark matter begins to be important. The dark matter annihi- 
lation rate per volume is proportional to n 2 , and hence the clustering of the dark matter 
roughly replaces this with the average of the density squared, (n x ), which can be much 
larger than the homogeneous part, n 2 The importance of this effect was recognized in 
Ref. [26J for the calculation of the diffuse gamma-rays from dark matter annihilation (see 
Refs. [6] for recent studies). In the context of cosmic reionization, the effect of dark 
matter clustering was pointed out in Refs. [131 [HJ [T6l [TT] . 

The effect of dark matter clustering is represented by the enhancement factor A 2 (z), 
given by [26J 




dMM 




dM 



(14) 
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Figure 3: (Top) TT spectrum of the CMB anisotropy. The black solid line shows the 
WMAP5 best fit curve without dark matter annihilation. Also shown are the cases of 
dark matter with annihilation cross section (crv) e + e - = 10~ 24 and 5 x 10~ 24 cm 3 s _1 with 
mass of m x =1 TeV. (Middle) The same, but for TE spectrum. (Bottom) The same, but 
for EE spectrum. 
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Figure 4: Constraint on the dark matter annihilation cross section with 99% CL. into 
e + e~ (dashed) and p + ji" (solid) as a function of the dark matter mass. 



where 



A 2 M (z) 



A v - Ir (z) h(x m )xl 



3 if^Xjji) 

Here dn(z)/dM is the comoving number density of collapsed objects with mass M, which 
can be evaluated by the Press-Schechter theory [27J[28]. A vir (z) denotes the virial over- 
density which, in the ACDM Universe, is given by [29] 



A vir (2;) 



18tt 2 + 82y - 39y 2 



(16) 



where Q m (z) = fi m (l + z) 3 /[Qa + ^m(l + z ) 3 } an d y = Q m (z 

P 

Rvir/a with i? v ; r being the virial 



1. I n (x m ) is defined as 

In{%m) — J m dxx 2 d n (x) with d(x) denoting the density profile of a dark mater halo, as 



p(r) = pd{r/a) with typical halo size a. Finally, x 
radius of the objects with mass M : 



47T 



Rl iT A viT (z)p m (z) = M. 



(17) 



For a given halo with mass M, the typical halo radius a may have some correlation with 
the virial radius R Y i T (M). We adopt a model of Ref. [30J. 

It is also noticed that the mass integral in Eq. (fT4j) has a lower cut (M cut ) below 
which the dark matter does not cluster, depending on the kinetic decoupling temperature 
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of the dark matter particle [211 E2], and this choice determines the redshift at which 
the clustering effects become important, as explicitly shown in Ref. [5]. We assume 
M cut = 1O~ 1O M in the following. Hereafter we adopt the Moore density profile [33] . 
where d(x) = x~ L5 (l + ar* 5 ) -1 . It is noticed that the clustering effect is more significant 
for halo density profile with steep core, such as the Moore's one. In order to include the 
enhancement effect, we only have to replace n^(z') with n 2 (z')[l + A 2 (z')} in Eqs. ffTTl) 
and (1T21) . 

Fig. [5] shows an evolution of the ionization fraction when no effects of dark matter 
annihilation is taken into account (black solid), when dark matter with (av) = 5 x 
10 -24 cm 3 s _1 and m x = 1 TeV are included (green dashed), and when the effects of dark 
matter clustering is also included (red dotted). In the last case, we have multiplied A 2 (z) 
by a factor 10 in order to clearly see the effects of dark matter clustering. In this figure 
we have assumed that the main annihilation mode is into e + e~. It is clearly seen that 
the dark matter clustering at low redshift enhances the annihilation rate and increases 
the ionization fraction. This effect, however, is found to be subdominant for the CMB 
anisotropy, since the main effect on the CMB comes from the energy deposition around 
the recombination epoch (z ~ 1000). A situation in which the ionization fraction is largely 
affected by the clustering effect at the low redshift is excluded even if the clustering effect is 
ignored. Therefore, for the purpose of constraining dark matter annihilation models from 
the observation of CMB anisotropy, taking only the homogeneous dark matter component 
is sufficient, in agreement with the statement in Ref. [T6] . 

5 Conclusions and discussion 

We have studied the effects of dark matter annihilation on the CMB anisotropy, motivated 
by the fact that recent measurements of anomalous comic-ray positron/electron fluxes 
can be explained by the contribution from dark matter annihilation with fairly large 
annihilation cross section. In contrast to many preceding works, we have taken into 
account all the relevant energy loss processes of electrons and photons and estimated how 
the injected energy goes into ionization, excitation and heating of atoms. 

This method is applied to the calculation of CMB anisotropy, and it is found that 
WMAP5 results give stringent bound on the annihilation cross section. This is much 
stronger than that from BBN [9] and other astrophysical bounds. Gamma-ray observa- 
tions also give stringent bound [HE], but it has rather large uncertainty due to the lack of 
the understandings of the dark matter halo profile. Taking in mind all the astrophysical 
uncertainties on these cosmic-ray signatures, the CMB constraint given in the present 
paper is said to be robust : the physics of CMB anisotropy is well understood from both 
theoretical and observational sides. Moreover, the forthcoming Planck satellite will give 
more precise data of the CMB anisotropy including polarization toward much smaller 
scales than the WMAP, and hence the bound is expected to be significantly improved. 
Furthermore, the heating of the intergalactic medium by the dark matter annihilation 
may have an impact on the 21cm observations [MJ [35]. We will pursue these issues in a 
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Figure 5: Evolution of the ionization fraction when no effects of dark matter annihilation 
is taken into account (black solid), when dark matter with (av) = 5 x 10~ 24 cm 3 s _1 and 
m x = 1 TeV are included (green dashed), and when the effects of dark matter clustering 
is also included (red dotted). We have assumed that the main annihilation mode is into 
e + e~. 



separate publication. 
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